The laser polarization as control parameter in the pattern formation 
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The recently observed dependence of the periodic surface structures on the light polarization 
in the laser induced pattern formation is analyzed within a model where the polarization induces 
significant deviation the spatial distribution of the energy deposited by the photon from isotropic 
energy distribution. We argue that the laser polarization breaks the rotation symmetry on the 
surface and is responsible for the correlation of the surface structures with the degree and the 
direction of polarization. Moreover it is shown that the polarization induces the appearence of novel 
features of the surface morphology and time evolution, which could be directly tested experimentally. 

PACS numbers: 74.25 



The phenomenon of ripple formation on the surfaces of 
eroded materials is of particular interests due to the large 
area of application in physical, chemical and material sci- 
ences. The development of periodically modulated struc- 
tures as result of sputtering, the removal of atoms from 
the surface of solids through the impact of energetic par- 
ticles (photons or ions), was discovered and studied ex- 
perimentally nearly three decades ago. The first widely 
accepted theoretical approach 0, Q describing the pe- 
riodic surface structures induced by laser radiation sug- 
gested that the ripples are the result of interference of in- 
coming laser beam with some form of a surface-scattered 
electromagnetic wave. In general, this theory was suc- 
cesseful in discription of uniformly distributed patterns 
with the periodicity dependent from the wavelength of 
laser radiation and from the angle of incidence. 

The development of ultra-short pulses laser technolo- 
gies and application of femtosecond laser radiation in 
surface sputtering [s| has shown that newly rediscov- 
ered laser induced periodic surface structures (LIPSS), 
with lateral periods a few times smaller than the wave- 
length of the incident light, can not been described in the 
framework of the conventional LIPSS theory Thus, it 
turns out that the ripple structure has non-trivial surface 
morphology sharing many similarities with aeolian sand 
dunes [5(, with periodicity independent from laser wave- 
length and the angle of incidence, but correlated with 
the local intensity of the laser beam that gives the strong 
support to the nonlinear self-organized mechanism of for- 
mation this structures. In this light the understanding of 
the puzzling dependence the ripple orientation from the 
laser polarisation [3, [6| where the resulting ripple orien- 
tation depends on the direction of the vector of electro- 
magnetic field, become very important and interesting 
problem. 

A revival of interest in the ion induced ripple forma- 
tion has been caused by successful theoretical prediction 
7] the ripple wavelength and orientation in agreement 
with experimental observations. Further development 
the nonlinear theory 0, [ij [13, EH allowed the proper 
description the time evolution of the surface morphology 



under ion-bombardment, ripple stabilization, wavelength 
dependence with ion energy or flux and production of dot 
structures as a function of bombardment conditions. 

In this Letter we develop a continuum theory of ero- 
sion by polarized laser radiation. We exploit connections 
with ripple formation by ion-beam sputtering and extend 
this model with inclusion of laser polarization, leading to 
polarization dependence of coefficients in nonlinear equa- 
tion of the Kuramoto-Sivashinsky type. Our results sug- 
gests that the laser polarization can be very important 
control parameter. 

Following to the Sigmund's theory of sputtering [l2| 
we consider the normal erosion velocity at the surface 
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where Vl'(r) is the local correction to the uniform flux due 
to variation of the local slopes, A is a material constant 
and E(r) is the average energy deposited at point O due 
to the scattering of the photon flux at P 
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e is the total energy carried by the photon flux and a, 
[3 and 7 are the Gaussian distribution widths along x',y' 
and z' axes, respectively, in the reference frame of the 
incoming beam. According to Sigmund's theory the en- 
ergy distribution widths a, /3, 7 along x',y',z' axis, re- 
spectively, scaled asa = /3>7^a - average depth of 
energy deposition. Obviously, the Gaussian form is not 
universal and the condition a = (3 is not a general case, in 
particular in the case of light interaction with a matter. 
Whereas the deviations from the Gaussian only slightly 
modified the pattern formation process (l3j , the asymme- 
try of energy distribution along surface plain has not been 
studied yet. One can expect, that the space anisotropy of 
the energy deposition can be effectively described by the 
anisotropy of the distribution widths a and (3. Therefore, 
we except the Gaussian form of energy deposition and 
consider the laser radiation which penetrates the bulk 
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of the material and stops at some point, where its en- 
ergy spread out to the neighboring sites (see Fig.l). In 
order to introduce the polarization in our calculations 
we choose the vector E of electromagnetic field to lie in 
the x'-z' plane in the reference frame of the incoming 
beam. Following to @, [1] we perform the calculation of 
erosion rate in the local coordinate system (X,Y,Z) shown 
in Fig.l. The local correction to the incident energy flux 
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FIG. 1: Reference frames for the computation of the erosion 
velocity: (x , y , z ) is reference frame of incoming beam where 
x' axis is parallel to the vector E of electromagnetic field; 
(X, Y, Z) correspond to the local coordinate frame, where 
Z axis is parallel to the local normal to the surface, while 
(x, y, h) denotes the laboratory frame of coordinates with h 
axis perpendicular to the flat substrate. The dotted ellipses 
are Gaussian distribution of deposited energy with halfwidth 
a,/3, 7 along x',y',z' axis, respectively. The incidence angle 
measured in the local reference frame is ip, and 9 in the lab- 
oratory frame. <j> is the angle between the incidence plane 
(X, Z) and E. Insert: view from top along h for 9 — 0. 

is now given by \&(r) = £Io(cosLp + (dxZ) siny>), where £ 
is the absorption coefficient of the material and Iq is the 
laser intensity. In order to describe the surface profile in 
the neighborhood of O we took into account cross-terms 
of the type ~ XY 
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As in [8[ we assume that the radii of curvature 
Rx,Ry,Rxy of the surface are much larger than the 
penetration depth a [l8j], so that only terms up to first 
order in a/ Rx, a/Ry and a/Rxy are kept. The integra- 
tion results in the erosion velocity V(ip, <j>, Rx , Ry , Rxy) 
as a function of angles tp, 4> an d the curvatures 1/ Rx — 
-d 2 Z/dX 2 , 1/Ry = -8 2 Z/dY 2 and 1/R XY = 
—d 2 Z/dXdY. Next we perform expansions in powers 
of derivatives of h(x,y,t) and rewrite V in terms of the 



laboratory coordinates (x,y,h) Q. 



dh(x,y,t) 
dt 



-v[$, <p, R x , Ry, RxyW 1 + (V/i) 2 (4) 



We complete Eq. (Q} by adding the following physical 
processes: the surface self-diffusion effects — JfV 2 (V 2 /i) 
where K is the relaxation rate due to thermally activated 
surface diffusion, together with the fluctuations r](x, y, z) 
(short noise) in the flux of the bombarding particles. Fi- 
nally, we obtain the equation of motion known as an 
anysotropic noisy Kuramoto-Sivashinsky [l^ | equation 
where the coefficients are now the complex functions of 
two angles (f> and 9. In order to simplify our consideration 
and because of the reasons discussed below we write here 
the equation for the case of the normal incident (6 = 0) 
and in the reference frame rotated by the angle <f> (see In- 
sert Fig.l) that means that v xy = [y x — u y ) tan(2<^>) = 
(X x - A y )tan(20) = 
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where the coefficients are given by 
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[b = a, /3, 7). Moreover, 



we neglected here the erosion velocity vq which does not 
affect the ripple characteristics, such as ripple wavelength 
and the ripple amplitude, and can be eliminated by the 
transformation h = h + Dot. 

The basic role of the introduced above parameters A 
and II can be clarified by means of the linear stability 
analysis of equation ([S]). The Fig. 2 shows the linear 
growth rate along q x and q y for various values of A and 
II. For A = the uniform state is stable, whereas for 
nonzero A and II = the instability sets both along q x 
and q y . An increase of II induces the asymmetry in q xy - 
plain thus for II ~ 1 the instability along q y is suppressed. 
Thus, the deviation from isotropic energy distribution 
breaks the rotation symmetry in xy plain and lids to an 
anysotropic linear instability. For the sake of definitcness 
we will fix A = 1/2 that corresponds to the condition 
a 2 + f3 2 — 7 2 . From ([5]) follows that the parameter II 
contributes to the both tension coefficients 
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One can see that the surface tension coefficients are nega- 
tive for the normal incidence and in general are not equal 
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FIG. 2: Stability of Fourier modes for equation ((Sj as a 
function of the wave vector q x (right) and q y (left), for various 
values of the parameters:A = 0,11 = (solid line); A = 1,11 = 
(dashed line); A = 1,11 = 0.4 (dotted line); A = 1,11 = 0.9 
(dash-dotted line). Insert: l x /l y as a function of II. 

to each other due to the fact that the direction of the 
laser polarization breaks the symmetry along the surface. 
Consequently, we expect the instability to the ripple for- 
mation with wavelength U = 2ny / 2K/\h'i\, where i refers 
to the direction (x or y) along which the associated [y x 
or v y ) is largest. Thus, in our case of polarization along 
x-axe for v x < v y < which holds when < II < 1 the 
ripple structure is oriented in the x direction. Moreover, 
the ripple wavelength has now the following dependence 
from IT 

^-^Miin) (8) 

Thus, for n = the wavelength l x = l y . The increase of 
LT slightly reduces l x whereas l y become to be very large 
for n — > 1 (linear polarization along x axis). The inset 
of Fig.2 shows the relation l x /l y = - U)/(l + II). 
Thus, the parameter II and the angle <fi allow us to change 
the orientation and topology of ripples that means that 
the parameter II is a control parameter which we will call 
in the following the degree of polarization. We will con- 
sider three cases: (i) circular polarization, II = (a a = 
cr/3); (ii) linear laser polarization, II = 1 (a a << crp), and 
(iii) elliptic polarization, < II < 1 (<r a < crp). 

Another quantity which can change the orientation of 
ripple is incidence angle 9. Thus, for angles 9 less than 
a critical angle 9 C , the wave vector of the modulations 
is parallel to the component of the beam in the surface 
plane |7(. By this means the orientation of ripple due 
to polarization and due to incoming angle can compete 
or work together. Fig. 3 displays the ripple orientation 
phase diagram (8, II) for three values of parameter o 1 



FIG. 3: Ripple orientation phase diagram for two val- 
ues of the parameter <r 7 — 1, 1.5, 2 (dashed, dotted, solid 
line, respectively). Region I: v x < v y \ Region II: v x > v y . 
Inserts show the incidence angle, the direction of ripples and 
the vector E. 

(negative values of II correspond to the vector E along y 
axis). The boundary is defined by ^(11, 9) = Vy(Tl, ff) 
and separates the region I (y x < v y ) and region II 
(y x > Vy). In the case of II = (circular polarization) 
and the normal incidence of laser beam the dot struc- 
tures are expected in our model. The increase of II for 
small aligned the ripples along the vector E. Whereas 
the growth of 9 can still change the direction of ripple for 
II > at some 9 C , which depends from er 7 , the laser po- 
larization completely suppress the transition for negative 
II. The interesting result, which can be tested exper- 
imentally, is the possible reorientation of ripple in the 
case of weak polarization —0.4 < II < for er 7 = 2. In 
any case, one have to perform detailed experimental in- 
vestigation the angle dependence for various values of II 
in order to fix the parameter er 7 . Moreover, our model 
allows us to understand the nature of LIPPS which are 
oriented parallel to the laser p-polarization at a relatively 
large incident angle. Nevertherless, we have to note here 
that our approximation may be not valid for the large 
angles of the incidence 9 where the effects of interaction 
of laser radiation with the material surface can not be 
neglected. Therefore, we limit ourselves in the following 
to the presentation of our results only for the normal in- 
cidence of laser beam and set a 1 = 2. The dependence of 
ripples structures on the angle of the incidence for large 
9 will be discussed somewhere else. 

Having analyzed the linear regime we now precede by 
investigating the influence of nonlinear terms and X y . 
As shown early [2(|, there is a clear separation of the 
linear and the nonlinear behaviour in time in such way 
that in the linear regime up to a crossover time t c the 
nonlinear terms as would be absent whereas the nonlinear 
terms take over after t c ~ {K / v 2 )ln(v / X) and completely 
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FIG. 4: Grey-scale plot of a surface, showing the surface 
morphologies for (from top to bottom): (j> — 0,11 = 1; 
4> = 30°, II = 0.5; 4> — 45,11 = 0, and for various values 
(from left to right) of time: t « t c ; t = t c /2 and t ~ t c 



normal parallel to the laser polarization, we consider the 
energy deposition along the laser polarization due to the 
powerflow distribution an oblate nanoparticle [13] ■ Thus, 
the concentration of the powerflow along the vector E of 
electromagnetic field near the nano-particles leads to the 
local instability of the lattice that could bring a collapse 
of the atomic structure about. We think that the surface 
atoms of such annealed areas have a good chance to be 
sputtered at the time scale of electron-phonon relaxation. 

In summary, we have shown that the correlation of 
ripples orientation with laser polarization can be con- 
sistently described within a model where the polariza- 
tion induces an anisotropy in the energy distribution and 
causes the symmetry break on the surface. The model ac- 
counts for experimental features of laser induced surface 
modulations and leads to the numerous predictions which 
can be directly tested experimentally. Moreover, our re- 
sults support the nonlinear self-organized mechanism of 
formation the ripples on the surface of solids. However, 
a more detailed investigation of the effective mechanism 
which is responsible for the anisotropy of energy transfer 
remains an interesting issue for future work. 



determine the surface morphology which depends on the 
relative signs of and X y . For our choice a 1 = 2 both 
X x and X y are negative and depend from II, that means 
the disappearance of ripples above t c ~ t c (H — 0)/(l + 
II) 2 , which is decreased with increasing of II and the 
appearance of kinetic roughening. 

The surface morphologies obtained by numerically in- 
tegrating the equation (0 for various values of 0, II 
and time show in Fig. 4 that the anisotropy in the energy 
deposition can describe the correlation of ripple orien- 
tation with laser polarization and saturation of ripples 
with increasing of time. However, in spite of remarkable 
concurrence the calculated and the ablated surfaces @ , a 
direct comparison with experimental date is difficult be- 
cause the relevant mechanism which is responsible for the 
energy anisotropy is still under discussion. Whereas in 
metalls and semiconductors the incident laser radiation 
can excite the plasma wave in the electron plasma pro- 
duced by multiphoton ionization [14J, the electron den- 



sity in dielectrics is too low for the intensities near the 
damage threshold. Instead, it was proposed, that the 
femtosecond laser pulses could induce high inhomoge- 
neous ionization under laser radiation inside the trans- 
parent solids and the dielectrics, forming nano-droplets 
of plasma by means of the " forest-fire" multiphoton and 
avalanche ionization 15j. Interaction of infrared laser 



radiation with nano-plasma-dielectric composites can ex- 
cite the surface plasmons polaritons (SPP) in the nano- 
particles and induce the giant enhancement of local elec- 
tric fields. In contrast to jl(f], where the underdense 
nanoplasmas grows into nanosheets orientated with their 
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